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We investigate here quantum effects in gravitational collapse of a scalar field model which clas- 
sically leads to a naked singularity. We show that non-perturbative semi-classical modifications 
near the singularity, based on loop quantum gravity, give rise to a strong outward fiux of energy. 
This leads to the dissolution of the collapsing cloud before the singularity can form. Quantum 
gravitational effects thus censor naked singularities by avoiding their formation. Further, quantum 
gravity induced mass fiux has a distinct feature which may lead to a novel observable signature in 
astrophysical bursts. 
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PACS numbers: 04.20Dw,04.60.Pp 

Naked singularities are one of the most exotic objects 
predicted by classical general relativity. Unlike their 
black hole siblings, they can be in principle directly ob- 
served by an external observer. There have been many 
investigations which show that given the initial density 
and pressure profiles for a matter cloud, there are classes 
of collapse evolutions that lead to naked singularity for- 
mation (see e.g. for some recent reviews), subject to 
an energy condition and astrophysically reasonable equa- 
tions of state such as dust, perfect fluids and such others. 
This has led to extensive debates on their existence, with 
a popular idea being cosmic censorship conjectures which 
forbid classical nakedness ^j. Since naked singularities 
originate in the regime where classical general relativ- 
ity is expected to be replaced by quantum gravity, it 
has remained an outstanding problem whether a quan- 
tum theory of gravity resolves their formation. Also, with 
the lack of observable signatures from the Planck regime, 
naked singularities could in fact be a boon for a quantum 
theory of gravity. Because, the singularity being visible, 
any quantum gravitational signature originating in the 
ultra-high curvature regime near a classical singularity 
can in principle be observed, thus providing us a rare 
test for quantum gravity. 

One of the non-perturbative quantizations of gravity is 
loop quantum gravity 0] whose key predictions include 
Bekenstein- Hawking entropy formula 0. Its application 
to symmetry reduced mini-superspace quantization of ho- 
mogeneous spacetimes is called loop quantum cosmol- 
ogy 5] whose success includes resolution of the big bang 
singularity Q , initial conditions for inflation 0, Is] j and 
possible observable signatures in cosmic microwave back- 
ground radiation These techniques have also been 
applied to resolve black hole singularity in a scalar fleld 
collapse scenario 0- 

Since the dynamics of a generic collapse is very com- 
plicated and tools to address such a problem in quantum 
gravity are still under development, it is useful to work 
with a simple collapse scenario as of a scalar fleld. It 
serves as a good toy model to gain insights on the role 



of quantum gravity effects at the late stages of gravita- 
tional collapse. Existence of naked singularities in these 
models is well-known and one of the simplest setting 
is to consider an initial conflguration of a homogeneous 
and isotropic scalar field <I> = <I>(t) with a potential V{^) 
(given by eq.©) and the canonical momentum P$. In 
this case it has been shown that fate of the singularity be- 
ing naked or covered depends on the rate of gravitational 
collapse 01 . For an appropriately chosen potential, for- 
mation of trapped surfaces can be avoided even as the 
collapse progresses, resulting in a naked singularity with 
an outward energy flux, in principle observable. Since the 
interior of homogeneous scalar fleld collapse is described 
by a Friedmann-Robertson- Walker (FRW) metric, tech- 
niques of loop quantum cosmology can be used to inves- 
tigate the way quantum gravity modifles the collapse. 

Let us consider the classical collapse of a homogeneous 
scalar field $(t) with potential V{^) and the canonical 
momentum for the marginally bound {k — 0) case. The 
interior metric is given by 



\t) [dr 



^dn^ 



(1) 



with classical energy density and pressure of the scalar 
fleld, 

p(t)-$V2 + V^('i>), Pit) = ^y2-Vi^) . (2) 

The dynamical evolution of the system is obtained from 
the Einstein equations which yield 0| 

R^R = F{t, r), p = F,r/KaR^, p = -F/kR^R (3) 

Here k = SttG, and F{t,r) — {k / i) p{t)r^ has interpre- 
tation of the mass function of the collapsing cloud, with 
F > and R{t,r) = ra{t) is the area radius of a shell 
labeled by comoving coordinate r. In a continual collapse 
the area radius of a shell at a constant value of comoving 
radius r decreases monotonically. The spacetime region 
is trapped or otherwise, depending on the value of mass 
function. If F is greater (less) than R, the the region 



2 



is trapped (untrapped). The boundary of the trapped 
region is given hy F = R. 

The collapsing interior can be matched at some suit- 
able boundary r = rh to a generalized Vaidya exterior 
geometry, given as HjI , 



and shall be constrained by phenomenological consider- 
ations. 

The change in behavior of the classical geometrical 
density (1/a^) for scales a < a,, can be well approxi- 
mated by I3 



ds^ = -il-2M{r^,v)/r,)dv^-2dvdr^+rldn^ . (4) d,{a) = D{q) , q:=a''/al a, := /jtTs^p (7) 

The Israel-Darmois conditions then lead to 0, 0| 
rba{t) = ry{v), F{t,rb) = 2M{ry,v) and 



with 



M{ry,v),r,, = F/2ri,a + rlaa 



(5) 



The form of the potential that leads to a naked sin- 
gularity is determined as follows. The energy density 
of scalar field can be written in a generic form as p = 
l"^^a~", where n > and I is a proportionality con- 
stant. Using energy conservation equation, this leads to 



the pressure p = [{n — 3)/3] 
eq.(I21) in these we obtain [llj 



in— 4. 



On subsituting 



-4:„2-n-2 



(6) 



Thus 



V{<P) = (1-71/6)/" 

Then it is easily seen that F/R = (k/3)/" 
in the collapsing phase as a — > 0, whether or not the 
trapped surfaces form is determined by the value of n. 
It is straightforward to check that for < n < 2, if no 
trapped surfaces exist initially then no trapped surfaces 
would form till the epoch a{t) — with a{t) = 

(l-nV2V3)'^". 

The absence of trapped surfaces is accompanied by a 
negative pressure implying that for a constant value of 
the comoving coordinate r, F is negative and so the mass 
contained in the cloud of that radius keeps decreasing. 
This leads to a classical outward energy flux. As the col- 
lapse proceeds, the scale factor vanishes in finite time and 
physical densities blow up, leading to a naked singularity. 
Since no trapped surfaces form during collapse, the out- 
ward energy flux shall in principle be observable. How- 
ever, near the singularity when energy density is close to 
Planckian values, this classical picture has to be modi- 
fied and we need to investigate the scenario incorporating 
quantum gravity modifications to the classical dynamics. 

Let us hence consider the non-perturbative semi- 
classical modifications based on loop quantum gravity 
for the interior. The underlying geometry for the FRW 
spacetime in loop quantum cosmology is discrete and 
both the scale factor and the inverse scale factor op- 
erators have discrete eigenvalues 0. In particular, 
there exists a critical scale a, — jj/Mp below which 
the eigenvalues of the inverse scale factor become pro- 
portional to the positive powers of scale factor. Here 
7 « 0.2375 is the Barbero-Immirzi parameter 0|, £p is 
Planck length and j is a half-integer free parameter which 
arises because inverse scale factor operator is computed 
by tracing over SU(2) holonomies in an irreducible spin j 
representation. The value of this parameter is arbitrary 



(,) = (8/77)%-V2{7[(, + l)"/4„|,„i|ii/4 

{q + l)y'-sgn{q-l)\q-l\y^]y. (8) 
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For a <C a*, dj cc [a/ a^,) and for a ^ a* it behaves 

classically with dj k, a^^ . The scale at which transition 
in the behavior of the geometrical density takes place is 
determined by the parameter j . 

At the fundamental level the dynamics in the loop 
quantum regime is discrete, however, recent invsetiga- 
tions pertaining to the evolution of coherent states have 
shown that for scales ao = ^/^(-p ^0,^0.* — \/J7/3£p, 
dynamics can be described by modifications to Fried- 
mann dynamics on a continuous spacetime [l^ with the 
modified matter Hamiltonian 



=dj{a)Pi/2 + a^V{'S>) 
and the modified Friedmann equation 

d^/a^ = {k/3){^^/2D + V{<i>)) 



(9) 



(10) 



which is obtained by the vanishing of the total Hamil- 
tonian constraint and the Hamilton's equations: $ = 
dj{a)P^, P$ = -a3T/$(cf>) 7]. These also lead to the 
modified Klein-Gordon equation 



+ {ia/a- b{q)/D{q)) ^ + D{q)V^^{<^) = 0. (11) 



Since at classical scales (a » a,) £> « 1, the modified 
dynamical equations reduce to the standard Friedmann 
dynamical equations. For scales a < a*, the $ term acts 
like a frictional term for a collapsing phase. We note that 
since semi-classical modifications for inhomogeneous case 
are still not known, we cannot do a complete quantum 
analysis of interior and exterior. The exterior is assumed 
to remain classical. Further, as a continuous spacetime 
can be approximated till scale factor ao, the matching of 
interior and exterior spacetimes remains valid during the 
semi-classical evolution. 

The modified energy density and pressure of the scalar 
field in the semi-classical regime can be similarly obtained 
from the eigenvalues of density operator and using the 
stress-energy conservation equation [3| 



Peff = (a) = $ V2 + D{q) 



(12) 
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FIG. 1: Evolution of area radius with time. The classical 
evolution (dashed) leads to naked singularity in finite time 
whereas in semi-classical evolution (solid) it is avoided. Inset 
shows evolution of energy density (in Planck units) with time. 
The parameters chosen are n — 1.9 and j = 100. 



and 



Poff 



2 1 D{q) 

3 (a/a) D{q) 
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D[q)V{^)-^^Vm. 
6{a/a) 

(13) 

It is then straightforward to check that PcS is generically 
negative for a < a, and for a ^ a* it becomes very 
strong. For example, at a ~ ao, Pcff ~ ~9poff- This is 
much stronger than its classical counterpart p — [{n — 
3)/3] p with < n < 2. Thus we expect a strong burst of 
outward energy flux in the semi-classical regime. Further, 
for a ^ a*, D{q) ^ 1 and the Klein-Gordon equation 
yields <i> ex: a^^ . Hence from the eq. (|12() we easily see 
that the effective density, instead of blowing up, becomes 
extremely small and remains finite. 

The modified mass function of the collapsing cloud can 
be evaluated using eq.® and ca. lfTUI) . 



(14) 



In the regime a ^ gq, d- becomes proportional to 
a^^, the potential term becomes negligible and thus the 
mass function becomes vanishingly small at small scale 
factors. 

The picture emerging from loop quantum modifica- 
tions to collapse is thus following. 

• Before the area radius of the collapsing shell reaches 
i?* = ra* at i = , collapse proceeds as per classical dy- 
namics and as smaller scale factors are approached $ and 
the energy density p cx a~" increase. The mass function 
is proportional to a"~'^ and (as < n < 2) it decreases 
with decreasing scale factor so there is a mass loss to 
the exterior, which is also understood from existence of 
negative classical pressure. 

• As the collapsing cloud reaches i?*, the geometric 
density classically given by a~^, modifies to dj and the 



dynamics is governed by the modified Friedmann and 
Klein-Gordon equations. The scalar field which experi- 
enced anti-friction in classical regime, now experiences 
friction leading to decrease of <!>. 

• The slowing down of $ decreases the rate of collapse 
and formation of singularity is delayed. Eventually when 
scale factor becomes smaller than oq this leads to break- 
down of continuum spacetime approximation and semi- 
classical dynamics. Discrete quantum geometry emerges 
at this scale 0| and the dynamics can only be described 
by quantum difference equation. The naked singularity 
is thus avoided till the scale factor at which a continuous 
spacetime exists. 

We show the evolution of area radius in time as collapse 
proceeds in Fig^ The semi-classical evolution (solid 
curve) closely follows classical trajectory (dashed) till the 
time f*. Within a finite time after t*, the classical col- 
lapse leads to a vanishing R and naked singularity. How- 
ever, the area radius never vanishes in the loop modified 
semi-classical dynamics and the naked singularity does 
not form as long as the continuum spacetime approxi- 
mation holds. The inset of FigQ] shows the evolution of 
energy density in Planck units. Classical energy density 
(dashed curve) blows up whereas it remains finite and in 
fact decreases in the semi-classical regime. 

The phenomena of delay and avoidance of the naked 
singularity in continuous spacetime is accompanied by a 
burst of matter to the exterior. If the mass function at 
scales a ^ a* is Fi and its difference with mass of the 
cloud for a < a* is AF = Fi — F, then the mass loss can 
be computed as 



AF 

F{ai) 



Pcsd, 



I-' 



(15) 



For a < a*, as the scale factor decreases, the energy den- 
sity and mass in the interior decrease and the negative 
pressure strongly increases. This leads to a strong burst 
of matter. The absence of trapped surfaces enables the 
quantum gravity induced burst to propagate via the gen- 
eralized Vaidya exterior to an observer at infinity. The 
evolution of mass function is shown in Fig|21 In the semi- 
classical regime, AF/Fi approaches unity very rapidly. 
This feature is independent of the choice of parameter 
j. The choice of potential causes mass loss to exterior 
in classical collapse also, but it is much smaller and in 
any case the classical description cannot be trusted at en- 
ergy density greater than Planck, when we must consider 
quantum effects as above. 

Interestingly, for a given collapsing configuration, the 
scale at which the strong outward flux initiates depends 
on the loop parameter j which controls a*. If j is large 
then burst occurs at an earlier area radius and vice versa. 
The inset of Fig [3 shows the mass loss ratio for different 
values of j. For all choices, AF/Fi 1, but the out- 
going flux profile changes. The loop quantum burst has 
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FIG. 2: Evolution of mass function with area radius for same 
parameters as in FigQ Loop quantum evolution (solid) leads 
to dissolution of all the mass of the collapsing shell. Dashed 
curve shows classical trajectory. Inset shows mass loss profile 
for j = 10*^ (outer), j = 5.0 x 10^ (middle) and j = 10^ 
(inner) . 

a distinct signature, at a ~ a.^ the flux decreases for a 
short period and then rapidly increases. Since the causal 
structure of classical spacetime is such that trapped sur- 
face formation is avoided, this quantum gravitational sig- 
nature can be in principle observed by an external ob- 
server as a slight dimming and subsequent brightening of 
the collapsing star. This peculiar phenomena is directly 
related to the peak in the function dj(a), and depends 
solely on the value of parameter j. If we compare this to 
other phenomenological applications 

this effect 

could not be masked by the role of other loop quantum 
parameters in a more general setting. This phenomena 
is thus a direct probe to measure j and an observer can 
estimate the loop quantum parameter j by observing the 
flux proflle of the burst based on this mechanism and 
measuring the variation in luminosity of the collapsing 
cloud. 

During such a burst most of the mass is ejected and 
this may dissolve the singularity. Thus non-perturbative 
semi-classical modifications may not allow formation of 
naked singularity as the collapsing cloud evaporates away 
due to super-negative pressures in the late regime. It 
has been demonstrated that these super-negative pres- 
sures would exist for arbitrary matter configurations 
which implies that results obtained here would hold even 
in a more general setting 16]. Loop quantum effects 
then imply a quantum gravitational cosmic censorship, 
alleviating the naked singularity problem. We note that 
the semi-classical effects do not show that the singular- 
ity is absent, it is only avoided till scale factor oq, be- 
low which the semi-classical dynamics and matching may 
break down. If for a given choice of initial data, semi- 
classical dynamics is unable to completely dissolve the 



singularity, the final fate of naked singularity must be 
decided by using full quantum evolution. Even in such 
cases we have valuable insights from semi-classical loop 
quantum effects with the possibility of phenomenologi- 
cally constraining the j parameter. 

In the toy model considered, we showed that the clas- 
sical outcome and evolution of collapse is radically al- 
tered by the non-perturbative modifications to the dy- 
namics. Our considerations are of course within the 
mini-superspace setting, and the general case of inhomo- 
geneities and anisotropics remains open. However, the 
possibility of such observable signatures in astrophysical 
bursts, as originating from quantum gravity regime near 
singularity is intriguing, indicating that gravitational col- 
lapse scenario can be used as probes to test quantum 
gravity models. 
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